Entropic enhancement of spatial correlations in a laser-driven Rydberg gas 
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In a laser-driven Rydberg gas the strong interaction between atoms excited to Rydberg states 
results in the formation of collective excitations. Atoms within a so-called blockade volume share a 
single Rydberg excitation, which is dynamically created and annihilated. For sufficiently long times 
this driven system approaches a steady state, which lends its properties from a maximum entropy 
state of a Tonks gas. Using this connection we show that spatial correlations between Rydberg 
atoms are controlled by the number of atoms contained within a blockade volume. For a small 
number the system favors a disordered arrangement of Rydberg atoms, whereas in the opposite 
limit Rydberg atoms tend to arrange in an increasingly ordered configuration. We argue that this 
is an entropic effect which is observable in current experiments. 

PACS numbers: 32.80.Eo, 61.20.Ne, 34.20.Cf 



Rydberg gases are currently intensely investigated as 
they constitute a versatile platform for probing strongly 
correlated phenomena in quantum many-body systems. 
The hallmark of many-body Rydberg physics is the ex- 
citation blockade, which inhibits the simultaneous laser 
excitation of two atoms within a given distance due to 
the strong interaction between Rydberg atoms [T]. This 
blockade effect and the emerging exclusion volume crit- 
ically determine the excitation dynamics of an ultracold 
Rydberg gas manifesting themselves in a strong suppres- 
sion of excitation probabilities |3] and non-Poissonian 
counting statistics H]. 

Since the Rydberg blockade defines an exclusion vol- 
ume around each excited atom it seems instructive to 
compare the situation in ultracold Rydberg gases with 
that encountered in classical hard objects models which 
arc often viewed as paradigmatic systems describing sim- 
ple liquids. In these systems the binary interaction po- 
tential is modeled by an impenetrable, finite-size core 
associated with each particle. The thermodynamic prop- 
erties of these simple exclusion volume models are dom- 
inated by entropy. Depending on the density of hard 
objects these systems may, in particular, exhibit density- 
density correlations over a considerable range or even un- 
dergo crystallization in higher spatial dimensions (see [S] 
and references therein). 

Yet, the Rydberg blockade differs in two main aspects 
from the situation found in hard objects models. First, 
the effect is coherent involving superpositions of many 
quantum states, i.e., all atoms contained within the ex- 
clusion volume share a single delocalized Rydberg exci- 
tation. Second, the blockade is a dynamical effect that is 
induced by the interplay between the resonant laser field 
and strong Rydberg-Rydberg interactions. Thus, exci- 
tations are continuously created and annihilated [3 [7] 
leading to the formation of a gas in which the number 
of Rydberg atoms is not conserved. This point appears 
to be critical for the emergence of long-range density- 
density correlations, as ordered phases in Rydberg gases 



have so far only been predicted for situations where states 
containing a fixed number of Rydberg atoms are deter- 
ministically excited [5UlO). For this dynamical crystal- 
lization technique a tailored sequence of laser pulses is 
used to connect an initial state without Rydberg excita- 
tions to a state with an ordered arrangement of Rydberg 
atoms. However, spatial correlations that extend well be- 
yond the blockade radius and emerge spontaneously (i.e. 
without specifically addressing the desired many-particle 
state) have, to the best of our knowledge, not been re- 
ported yet [TTHn]. 

In this work we address the question whether pro- 
nounced density-density correlations can spontaneously 
develop in strongly interacting Rydberg gases by analyz- 
ing their dynamics at long times. It has been demon- 
strated experimentally [3 [TH [13] and theoretically [TBl - 
[TO] that observables like the number of excited atoms 
reach a steady state after some time. Moreover, recent 
theoretical work on Rydberg lattice gas models has indi- 
cated that properties of the steady state can be under- 
stood from the perspective of thermodynamic ensembles 
[T^ . Here, we use this insight and argue that the equilib- 
rium state of a Rydberg gas shares important properties 
with a thermal maximum entropy state of a Tonks gas 
|20) . Connecting to this result and using purely entropic 
arguments we identify the figure of merit that controls 
the Rydberg density as the number of atoms contained 
within the blockade volume. The larger this quantity, 
the higher the density of excited atoms and the more 
enhanced the spatial correlations of the Rydberg atoms 
become. 

We demonstrate this for a one-dimensional system. 
This has the advantage that a number of analytical re- 
sults can be derived and numerically exact simulations of 
the quantum dynamics of 50 — 100 particles are feasible. 
We consider a gas with atomic line density a and large 
but fixed length L. For the sake of simplicity (the phys- 
ical mechanism discussed here does not depend on that) 
we consider atoms that are arranged on a regular lattice 



with spacing so that the number of lattice sites is 
aL. The density of Rydberg atoms is denoted as p. Due 
to the blockade the maximum possible Rydberg density 
is limited to Pmax ~ where /b is the blockade radius 
(cf. Fig. [T]), which we will also assume to be fixed in the 
following. The internal structure of the atoms is mod- 
eled by two states: ||) is the ground state and |t) is the 
Rydberg state. Both states are resonantly coupled by a 
laser of Rabi frequency J7. The interaction between two 
excited atoms separated by a distance r is given by the 
van-der-Waals potential V{r) = C^/r^ with dispersion 
coefficient Cg. The Hamiltonian of the system is 
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with Uk = + 1)/2 and the Pauli matrices ctq. For suf- 
ficiently large interaction the simultaneous excitation of 
nearby Rydberg atoms is strongly suppressed. The cor- 
responding exclusion or blockade radius 1^ is calculated 
by setting the collective laser coupling of atoms within a 
blockade radius fl\/a Ih equal to the interaction energy 
V{r) of two Rydberg atoms at a distance Zb TS" yielding 
«b = [C6/(l^v^)]2/". This length is in fact not sharp and 
there is a finite though small probability that two atoms 
at a distance closer than Zb are simultaneously excited. 
However, in experiments and in theoretical calculations it 
has been shown that the assumption of a sharp blockade 
radius accurately captures the physics of in an interact- 
ing Rydberg gas |l6j. Furthermore, we also neglect all 
residual interactions beyond Ih due to the tail of the van- 
der-Waals potential as they are not essential to obtain a 
qualitative picture. Later on, we will further elucidate 
the consequences of this spatial cut-off. Within these ap- 
proximations the Hamiltonian can be transformed to the 
form [HllTllSg 
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where the exclusion of Rydberg atoms is made manifest 
by the projection operator formed by the term inside the 
square brackets. This projector probes whether there 
are excitations within the exclusion length from a given 
atom. If so, it yields zero while it is one in the opposite 
case. 

The Hilbert space of the system consists of a number 
of subspaces which are not coupled by Hamiltonian 
We are interested here in the experimentally relevant sub- 
space which contains the configuration |0) = |m ... |), 
i.e. the state without any Rydberg atom, and all config- 
urations in which there are no two excitations within a 
blockade volume (orange ellipses in Fig.[T]). All these clas- 
sical configurations span the portion "Hdyn of the Hilbert 
space in which the system evolves. In order to count the 
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FIG. 1. Excitation blockade and Tonks gas. Ground state 
atoms within a blockade radius Ih from a Rydberg atom (red 
solid circle) cannot be excited to a Rydberg state. This defines 
a blockade volume of size 2Zb around each excited atom (or- 
ange ellipses), in which a single Rydberg excitation is shared 
between all particles. Certain sites may be located in the 
intersection of two blockade volumes (grey circles). In or- 
der to count the number of possible arrangements of Rydberg 
atoms it is convenient to map the system to a Tonks gas of 
non-overlapping hard rods (blue rectangles) of length 



number of these configurations - which equals the dimen- 
sion of Hdyn - it is convenient to map the physical situ- 
ation, where blockade volumes may overlap to a system 
of non-overlapping hard rods (blue rectangles in Fig.[T]). 
These hard rods have a diameter l\, and occupy al^ + 1 
lattice sites. Counting the number of ways for placing 
these hard rods reduces to the combinatorics of a lattice 
Tonks gas [23]. This immediately yields the number of 
configurations containing N Rydberg atoms 
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and the Hilbert space dimension dim'H, 

Y^(TL/(CT!b + l) ^ (i.T\ 
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In one and two dimensions it has been shown that 
Hamiltonians of the form ^ induce quantum dynam- 
ics that lead to a steady state in which all configura- 
tions spanning "Hdyn are populated with equal probabil- 
ity [T8l[T9]. The steady state is, therefore, given by a a 
microcanonical maximum entropy state in the sense that 
expectation values of observables can be calculated from 
a microcanonical partition function. The partition sum 
for the system under study here is in fact given by Eq. 
([3]). This is the connection of the long-time dynamics 
of a strongly interacting Rydberg gas with the thermal 
state of classical hard objects models. 

For our one-dimensional setup we can use this connec- 
tion to compute the steady state properties analytically. 
The number of Rydberg atoms (N) in equilibrium can 
be obtained from the distribution function Si(iV) using 
(N) = (dimHdyn)-^ j:"N^t'"^'^ N^l{N). In the large 
volume limit {L ^ Zb) this function is strongly peaked. 
Thus, the mean number of Rydberg excitations in the 
steady state is given by the position of this peak which 
can be calculated using 9ArlogSi(A^)| 
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suming further that the number of atoms within each 
blockade volume is large {alh ^ 1) we obtain the rela- 
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FIG. 2. Probability for finding configurations witii a given 
density p of Rydberg atoms for different values of cr Zb in a 
system of length L — 200 x Ih- For a small number of ground 
state atoms within a blockade radius (crib) the function peaks 
at small densities p <^ l^^ (see top left panel for three typical 
configurations). When alh grows the relative weight of config- 
urations with a large Rydberg density p ~ l^^ increases dra- 
matically. Here a huge number of configurations with nearly 
crystalline arrangement exist (see top right panel for three ex- 
amples). Hence, when all configurations are populated with 
equal probability, a state with high density and enhanced spa- 
tial correlations is entropically favored. 
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with the line density of Rydberg atoms p — {N)/L. 
ing the Lambert VF-function defined by W{x) e^'y^^ 
this can explicitly be computed: 

^ 1 W{ah) 
^ h\ + W{ahy 

An analysis of this expression shows that for fixed length 
L the density of Rydberg atoms in the steady state grows 
as the total atomic density and, therefore, the number al\, 
of atoms per blockade radius is increased. This behavior 
is shown in Fig. [2] where the graph depicts the proba- 
bility for finding configurations with a given density p of 
Rydberg atoms for different values of cr Zb in a system of 
large but fixed size L and fixed Ih L. 

In addition to the enhancement of the Rydberg density, 
the number fluctuations relative to the mean number of 
excited atoms decrease with increasing number of par- 
ticles in the blockade volume. To characterize this we 
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determine the Mandcl Q-parameter: Q 
In the vicinity of its maximum the function ([s]) can be ap- 
proximated by a Gaussian. The fluctuations in the num- 
ber of Rydberg atoms are then given by {N'^) — {N)'^ 
-1/ {d%logEL{N)\N={N)), yielding 



Q = [ply^f - 2p/b - [1 
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The Q-parameter decreases with increasing crZb, he. the 
distribution function depicted in Fig. [2] becomes increas- 
ingly sub-Poissonian. Since the blockade radius is held 




FIG. 3. Density-density correlation function obtained from a 
Monte-Carlo simulation for a sample of length L — 10 x Zb (pe- 
riodic boundary conditions) and atomic densities correspond- 
ing to alh = {f , 5, 30}. The larger the number of ground state 
atoms per blockade radius uh the more pronounced become 
the oscillations of the spatial correlations. 



constant and the density of Rydberg atoms increases this 
behavior of the fluctuations indicates an enhancement of 
spatial correlations between the excited atoms. 

The fact that pronounced spatial correlations emerge 
when just the total atomic density a is increased, is a di- 
rect result of the counting of all configurations which led 
to Eq. This is a somewhat counterintuitive finding. 
One might expect that the maximum of the entropy func- 
tion (cx log Si (A'')) lies at low Rydberg densities irrespec- 
tive of the lattice spacing, since there is a vast number of 
possibilities for exciting a relatively small number of Ry- 
dberg atoms on the lattice. However, it is actually also 
possible to gain entropy with nearly ordered configura- 
tions. To illustrate this, suppose we have a configuration 
in which n Rydberg atoms are already present. Then, the 
number of possibilities to place an additional excitation 
is given by the number of sites that are not enclosed by 
any blockade volume. Hence, for closest packing of the 
n Rydberg atoms the number of available sites is max- 
imized. This is because (n — l)alh blockaded sites are 
shared among the already excited atoms (grey sites in 
Fig.[l]showing the case of non-maximum overlap of block- 
ade volumes). Clearly, this dense packing of n Rydberg 
atoms increases the entropy for placing another one. This 
effect becomes more pronounced the larger the number of 
atoms contained in a blockade volume. Eventually this 
outweighs the tendency of the maximum entropy state 
to be formed by low density configurations. Thus, in 
summary, for small values of alh the steady state mostly 
contains low density configurations (cf. top left panel in 
Fig. [2]) while in the opposite limit large Rydberg densi- 
ties are favored (cf. top right panel in Fig. [2]) . The latter 
implies enhanced spatial correlations. 

To characterize these correlations in more detail, 
let us study the Rydberg-Rydberg correlation function 
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FIG. 4. Density-density correlation function for a sample of 
length L = 10 X (periodic boundary conditions) and an 
atomic density of aZb = 7. The crosses are obtained from 
a time-averaged quantum calculation using Hamiltonian ([2| 
and the empty state |0) as an initial state. The solid curve 
interpolates data obtained from a (classical) Monte-Carlo sim- 
ulation for the same parameters. 

3(2) (|x|) = (iV(0)7V(a;))/(7V)2. To calculate g^^Hkl) we 
use a classical Monte- Carlo method that generates a rep- 
resentative sample of classical configurations. In Fig. [3] 
we show data for a sample of length L = 10 x l^, with 
periodic boundaries and for various values of the atomic 
line density: cr^b = {1,5,30}. Note that because of the 
underlying lattice, x is only defined at positions where 
a X assumes integer values. We find that for an increas- 
ing number of ground state atoms per blockade radius 
the oscillations in the correlation function become more 
pronounced. This is consistent with the previously dis- 
cussed decrease of the Q-parameter (|6| . For alh ~> 1 the 
discreteness of the ground state gas can be neglected and 
the correlation function is well approximated by that of 
a Tonks gas [53] 

_ e-^(-'^)N/'^ ^ g(.,)(a?bXfc)'- 
Pk t[ 

with Xk = |a;|/^b — k and the unit step function 9{x). 

From this we see that the first peak decays exponentially 

with decay constant [^b W^(o'^b)]^^ and the height of the 
(2) 

peak is given by ffTo„ks(^b) = W{alh) + 1. This agrees 
with our simulation results. 

Let us finally show that the steady state of the quan- 
tum evolution is indeed well captured by the statistical 
treatment and that we can indeed expect entropy induced 
spatial correlations to emerge in the coherently evolving 
quantum system. We have performed a numerical propa- 
gation of the initial state |0) under the Hamiltonian ([2]). 
Due to the exponential scaling of the Hilbert space di- 
mension this can only be achieved within a rather lim- 
ited parameter range. We have chosen cr/b = 7, a system 
size of 10 X /b and a final propagation time of fit — 150. 
Here the system is well within the steady state regime. In 



Fig. [4] we show the density-density correlation function 
averaged over the time-interval 130 < fit < 150. The 
averaging removes small temporal fluctuations that are 
due to the finite size of the system (see Refs. [T7HIS])- 
The agreement of the data (crosses) to the results from 
the classical Monte-Carlo simulation (solid line) is good. 
The first two peaks are almost perfectly reproduced, with 
the first one showing the expected exponential decay and 
the correct height. The agreement becomes less accurate 
with increasing |a;| but pronounced oscillation in the spa- 
tial correlations are nevertheless visible also here. The 
discrepancy with respect to the Monte-Carlo result is due 
to finite size effects and the particular choice of the initial 
state as discussed in Refs. [ITl fT9]. 

One-dimensional Rydberg gases as discussed here can 
be realized in elongated atomic clouds with a transverse 
extension much smaller than ?b, see e.g. [15l[25]. Since Ih 
is typically on the order of 5 ... 10 /xm, currently achiev- 
able densities would allow values of trZb = 30 or even 
larger. Yet, as compared to an experimental realization, 
the simulations presented here were undertaken for ide- 
alized conditions, that is, (i) an equidistant uniform dis- 
tribution of atoms and (ii) no influence of the tail of the 
interaction potential. Both points can, in principle, be re- 
alized with available experimental technology. The first 
condition can be approximately achieved in an optical 
lattice [2S] , the second by using microwave dressing tech- 
niques to "chop off" the tail of the interaction potential 
[26| [27] . But, even for experiments in a continuous gas 
and without tailoring the Rydberg- Rydberg potential our 
analytical results remain qualitatively valid. 

In the disordered gas the atomic positions will vary 
from shot to shot, i.e., the number of atoms per blockade 
volume fluctuates which will lead to a higher Q parame- 
ter than predicted by Eq. (|6|. The significance of these 
fiuctuations, however, decreases with increasing atomic 
density a, so that for sufficiently large densities our lat- 
tice gas approach becomes more and more accurate. The 
tail of the van-der-Waals potential will start to play a 
role for Rydberg densities close to l//b- In this regime, 
the concept of a fixed blockade radius is no longer ac- 
curate, as atoms that are not enclosed by any "blockade 
volume" may still be blockaded by the addition of the 
potential tails of two nearby Rydberg atoms. This will 
lead to a saturation the maximum achievable enhance- 
ment of Rydberg- Rydberg correlations [IS]. Due to the 
short-range character of the van-der-Waals potential this 
saturation effect will, however, happen at comparatively 
large Rydberg densities. Most importantly, the inter- 
esting effect of entropically enhanced spatial correlations 
of Rydberg atoms is already significant at relatively low 
Rydberg densities and will, thus, be experimentally ob- 
servable. 

In one dimension correlations always decay exponen- 
tially. In the future it will therefore be interesting to 
study the discussed system in higher dimensions as here 
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long range order can emerge. This is difficult to treat 
theoretically but could be observable with current ex- 
periments together with appropriate imaging techniques 

[HiMlIsn]. 

Acknowledgements — We acknowledge fruitful dis- 
cussions with J. Evers, M. Garttner and J. P. Gar- 
rahan. This work was funded in part by EPSRC 
Grant no. EP/I017828/1 and Leverhulme Trust grant 
no. F/00114/BG. C.A. acknowledges support through a 
Feodor-Lynen Fellowship of the Alexander von Humboldt 
Foundation. 



[2 



13: 

[4; 

[6 

[7; 

[9 
[10 

[11 



M. D. Lukin, M. Fleischhauer, R. Cote, L. M. Duan, 
D. Jaksch, J. I. Cirac, and P. ZoUer, Phys. Rev. Lett. 87, 
037901 (2001) 

D. Tong, S. M. Farooqi, J. Stanojevic, S. Krishnan, Y. P. 
Zhang, R. Cote, E. E. Eyler, and P. L. Gould, Phys. Rev. 



Lett.| 93, 063001 (2004) 



K. Singer, M. Reetz-Lamour, T. Amthor, L. C. Marcassa, 
and M. Weideniuller, |Phys. Rev. Lett.| 93, 163001 (2004) 
T. C. Liebisch, A. Reinhard, P. R. Berman, and 
G. Raithel, Phys. Rev. Lett. 1 95, 253002 (2005) 
S. Torquato and F. H. Stilhnger, Rev. Mod. Phys.| 82, 
2633 (2010) 

U. Raitzsch, V. Bendkowsky, R. Heidemann, B. Butscher, 
R. Low, and T. Pfau, Phys. Rev. Lett. 100, 013002 (2008) 
M. Reetz-Lamour, T. Amthor, J. Deiglmayr, and M. Wei- 
demiiller, Phys. Rev. Lett. 100, 253001 (2008) 
T. Pohl, E. Demler, and M. D. Lukin, Pliys. Rev. Lett.l 
104, 043002 (2010) 

J. Schachenmayer, L Lesanovsky, A. Micheli, and A. Da- 
ley, New J. Phys 13, 059503 (2011) 

R. M. W. van Bijnen, S. Smit, K. A. H. van Leeuwen, 

E. J. D. Vredenbregt, and S. J. J. M. F. Kokkelmans, J. 
Phys. B 44, 184008 (2011) 

F. Robicheaux and J. V. Hernandez, Phys. Rev. A 72, 
063403 (2005) 



[12] S. Wiister, J. Stanojevic, C. Ates, T. Pohl, P. Deuar, 
J. F. Corney and J. M. Rost, Phys. Rev. A, 81, 023406 
(2010) 

A. Schwarzkopf, R. E. Sapiro, and G. Raithel, |Phys. Rev. 



[13 



[14 



[15; 



[16 
[17 

M 



[19 

[20 
[21 
[22 

[23 



[24 
[25 

[26 

[27 

[28 
[29 



Lett. 107, 103001 (2011) 
R. Heidemann, U. Raitzsch, V. Bendkowsky, B. Butscher, 
R. Low, L. Santos, and T. Pfau, Phys. Rev. Lett. 99, 
163601 (2007) 

R. Low, H. Weimer, U. Krohn, R. Heidemann, V. Bend- 
kowsky, B. Butscher, H. P. Biichler, and T. Pfau, Phys.| 



Rev. A 80, 033422 (2009) 
H. Weimer, R. Low, T. Pfau, and H. P. Biichler, Phys. 
Rev. Lett. 101, 250601 (2008) 

B. Olmos, M. Miiller, and I. Lesanovsky, New Journal of 
Physics 12, 013024 (2010) 



L Lesanovsky, B. Olmos, and J. P. Garrahan, Phys. Rev. 
TetEj lOS, 100603 (2010) 

C. Ates, J. P. Garrahan, and I. Lesanovsky, Phys. Rev. 
Lett, (accepted); preprint, arXiv:1108.0270 (2011) 

L. Tonks, Phys. Rev. 50, 955 (1936) 

L Lesanovsky .Phys. R ev. Lett. 1 06, 025301 (2011) 

S. Ji, C. Ates, and L Lesanovsky, |Phys. Rev. Lett. 1 107, 

060406 (2011) 

D. Chowdhury and D. Stauffer, Principles of Equi- 
librium Statistical Mechanics (Whiley-VCH, Weinheim, 
Germany, 2000) 

Z. W. Salsburg, R. W. Zwanzig, and J. G. Kirkwood, J. 
Chem. Phys. 21, 1098 (1953) 

M. Viteau, M. G. Bason, J. Radogostowicz, N. Malossi, 
D. Ciampini, O. Morsch, and E. Arimondo, Phys. Rev.| 
Lett. 107, 060402 (2011) 

H. P. Biichler, E. Demler, M. Lukin, A. Micheh, 
N. Prokof'ev, G. Pupillo, and P. Zoller, ,Phys. Rev. Lett., 
98, 060404 (2007) 

A. Micheh, G. Pupillo, H. P. Biichler, and P. Zoller, Phys. 
Rev. A 76, 043604 (2007) 

M. Garttner and J. Evers, private communation 



B. Olmos, W. Li, S. Hofferberth, and L Lesanovsky, Phys. 



Rev A| 84, 041607 (2011) 
[30J G. Giinter 



M. Robert-de-Saint- Vincent, H. Schempp, 
C. S . Hofmann, S. Whitlock, and M. WeidemiiUer, Physl 
Rev. Lett. 108, 013002 (2012) 



